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Abstract. We present a regular cubic lattice solution to Einstein field equations 
that is exact at second order in a small parameter. We show that this solution 
is kinematically equivalent to the Friedmann-Lemaitre-Robertson- Walker (FLRW) 
solution with the same averaged energy density. This allows us to discuss the fitting 
problem in that framework: are observables along the past lightcone of observers 
equivalent to those in the analogue FLRW model obtained by smoothing spatially 
the distribution of matter? We find a criterion on the compacity of the objects that 
must be satisfied in order for the answer to this question to be positive and given 
by perturbative arguments. If this criterion is not met, the answer to this question 
must be addressed fully non perturbatively along the past lightcone, even though the 
spacetime geometry can be described perturbatively. 



1. What is the fitting problem? 

Cosmology is unique among physical sciences for a certain number of reasons. First, the 
Universe is given once and for all, and there is no possibility to compare it to another 
Universe. This can usually be overcome by supposing that the initial conditions for 
the Big Bang model must be generic in some reasonable sense, or that some mechanism 
(e.g. inflation) is responsible for making them generic. Second, we are the only available 
observers in the Universe; there might be other observers, but we do not have access 
to their observations. The only piece of information we have on the Universe comes 
from our past lightcone, and a few local (geological) measurements on our worldline. 
Therefore, in general, one cannot rely purely on observations to fully determine the 
nature and dynamics of the Universe: one has to introduce extra assumptions on the 
theory of gravitation, the geometry of the Universe on large scales and the physical 
nature of its matter content. In the present paper, we will suppose throughout that 
gravity obeys the laws of General Relativity, and we will concentrate on the other 
two points: the geometry of the Universe and its matter content. In the standard 
model of cosmology, it is assumed that 'on average', on sufficiently large scales, the 
distribution of matter in the Universe is well described by a set of perfect fluids whose 
energy densities and pressures are locally homogeneous and isotropic; this results in 
Friedmann-Lemaitre-Robertson- Walker (FLRW) Universes with spatial sections that are 
maximally symmetric, i.e. determined entirely by their constant Gaussian curvature. 
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This assumption is based on the observed almost isotropy of the Cosmic Microwave 
Background around us, together with the Copernican Principle, and is usually called the 
Cosmological Principle. It is clear that it is an extrapolation outside our past lightcone, 
since the notion of average implicitly present in this principle tells us something about 
the spatial distribution of matter, starting from its distribution along our past lightcone. 
In a nutshell, the fitting problem [1] can be summarized by the question: does the 
effective FLRW model obtained by extrapolating the observed properties down our past 
lightcone coincide with the effective FLRW model obtained by smoothing the spatial 
distribution of matter? Of course, this question is not independent on the set of observers 
used to define the notion of spatial distribution: it makes use of a preferred set of 
observers, called fundamental in the standard model; usually, in the late-time Universe, 
the fundamental observers are supposed to be comoving with virialised objects such 
as galaxies, so as to include us among fundamental observers. In this paper, we will 
try and address the fitting problem by considering a special dynamical solution to the 
field equations consisting in a regular cubic lattice of initial cell size L and of objects of 
equal masses M. The solution is exact at second order in the small parameter MjL 
and we will see that it exhibits, on average, the same dynamical behaviour than the 
equivalent flat FLRW model with a non-relativistic fluid of density p = M/L^, therefore 
showing that, at second order in ^J1<A~[L^ there is no backreaction in the model [2]. Then, 
because the solution for the metric is exact at this order, we will be able to solve the 
Sachs equations at second order in MjL in order to reconstruct observables such that 
the distance- redshift relation. We will see that this solution for observables presents 
some divergences linked with the compacity of the object: if the extension ri of the 
objects is too small, the perturbative expansion of the solution of the Sachs equations is 
no longer valid, even though the perturbative expansion of the solution of the Einstein 
field equations remains stable. Namely, we will show that observables in this model 
remain very close to the observables calculated in the analogue FLRW model with energy 
density p = M/L^, as long as the parameters of the lattice obey: . If 

this condition is not satisfied, then observables cannot be calculated perturbatively, even 
though the metric is well approximated by the perturbative expansion, and one must 
solve the full system of Sachs equations [3]. Results regarding this complete integration 
are hinted at in this paper. These results illustrate the importance of the fitting problem 
in cosmology: the kinematically averaged FLRW model and the FLRW reconstructed 
by fitting observations might differ significantly (if calculated perturbatively, at least), 
even in the absence of (kinematical) backreaction. We work in units G = c = 1. 

2. A lattice Universe: kinematics and observables 
2.1. The cubic lattice solution 

Let us start by describing the lattice solution (see [2] for a complete derivation); we will 
only sketch the results and discuss their implications. We start with a cubic lattice of 
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size L with identical masses M at the centre of each cell. If the masses on the lattice 
are to represent typical galaxies, we can choose, as our typical parameters M ~ lO^^M© 
and L ~ 1 Mpc, where Mq ~ 10^° kg is the Solar mass, and L is of order of the 
intergalactic distances. Then, the natural parameter of the lattice is Rs/L ~ 10~^ ^ 1, 
where Rs = 2M is the Schwarzschild radius of the masses. Therefore, we can look for a 
solution expanded into powers of y/W/L (in |2] we prove that there is no perturbative 
solution in powers of M/L); this will lead to linearised field equations that can be 
solved exactly. We choose coordinates that are comoving with the masses: 5^00 = — 1, 
and spatial coordinates are Cartesian and adapted to the symmetries of the lattice. 
The distribution of matter is therefore a three dimensional Dirac comb with the masses 
located at Xn = Ln, n e Z^; the energy momentum tensor is then: Tab = TqqS^S^ such 
that: 

Too ocM E (5(^)(x-Ln). 

Actually, the field equations do not have a solution for such a source term [H |2] , because 
the formal series solution presents a UV divergence coming from the point-like nature of 
the masses. Therefore, we introduce a UV cut-off by giving a small but finite extension 
to the masses, rj and by replacing the Dirac deltas by their standard approximation: 

6{x — uL) ~ • Then, we write the source term in Fourier series, and we 

expand the field equations in powers of a/ M/ L and solve them order by order, to find 
the following solution at second order: \/i G {1,2,3}, (^oi = 0, and W{i,j) G {1,2,3}^: 



9ij = 



^ ^ M 8n t 2M / , , 2nt^ 



where e = ±1 and: 



+ ^mu^) (1) 



Let us insist on the fact that this solution is exact at order M/ L. We can now calculate 
the rate of expansion between two masses of the lattice. For that, consider two masses 
on the X-axis (all the other axes are equivalent, by symmetry), separated by a coordinate 
distance NL, for an integer. The physical distance between the two masses is given by 
l{t) = j^^ y/Qxxdx, and, expanding the square root to order M/L we find the effective 
scale factor of the lattice: 

a(t) = — = 1 + e\ —\ —t — . (3) 

^ ^ A^L V 3 V 3L3 ^ ^ 

The Hubble flow defined by H{t) = a{t)/a{t) is then found to be, at order M/L: 

H{t) = s-^f - Therefore, defining the initial Hubble rate = 

and choosing the expanding solution, e = 1, we get: 

H{t) = Ho- iHlt + 0{Hl) (4) 
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and this corresponds exactly, at order M/L oc Hq to a flat FLRW model fllled with 
non-relativistic dust. The result is actually valid at order (M/L)^^^ [2J. Thus, the 
model with discrete masses on a cubic lattice, once smoothed, is identical to a FLRW 
model with dust, with the corresponding smeared energy density. This means that, 
from purely kinematical considerations, one cannot distinguish between the average, 
homogeneous fluid description of the lattice and the exact behaviour of this lattice: 
there is no backreaction (in the sense of [5]) associated with spatially smoothing the 
lattice. 



2.2. Ohservahles, compacity and the fitting problem 



Now that we have a solution of the fleld equations that does not display backreaction, 
we can try and address the fltting problem by comparing observables in the lattice with 
observables in the corresponding, smoothed FLRW Universe with the same kinematics: 
any discrepancy between the two will be a sign that there exists a fltting problem. In 
order to carry the comparison, we calculate the distance/redshift relation in the lattice. 
The 4- velocity of an observer in one of the objects of mass M is given by u'^ = (1, 0, 0, 0) 
according to our choice of coordinates, and we deflne A, an affine parameter down light 
rays. We denote by O and 5* the locations of observer and the source respectively. Given 
the normalisation chosen in |3], the null vector of a past-directed light ray, is such 
that kQ = 1, so that we have, for the redshift: 1 + z{X) = l^a"")'^ = (A). Therefore, 
we can solve the component of the null geodesic equation order by order in terms 
of a/M/L. The distance is obtained similarly by solving Sachs equations [6] expanded 
at order M/L. The Sachs equations are actually exactly solvable in this perturbative 
scheme because the equations for the isotropic expansion and the shear decouple from 
each other at that order. Details of the calculations can be found in [3] . In terms of the 
past-directed affine parameter A < (and A = at the observer), we find that, at order 
M/L: 



ziX) 




+ 



M 1871 X 

t\Iyl 



(5) 



x + 



271 M X^ 



1+ ^~ 



2M ^ 



X 



1=4 

E 

1=1 



-X 



cos (^^^) L sin (^^^) 



TT 



(6) 



Here: Ui = {n,p,q), U2 = {n,-p,-q), U3 = {n,p,-q), U4 = {n,-p,q), and 
= {{n,p,q) G : 3/ G {1, 2, 3, 4}/u;.v = 0}. This means that the first sum is 
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over all the triplets that cancel one at least of the u^.v, whereas the second sum is over 
all the other triplets. These expressions coincide with their FLRW counterparts for the 
FLRW model obtained by smoothing the distribution of masses of the lattice, up to 
a priori small corrections proportional to M/L (the parts that are non-polynomial in 
A). Actually, it turns out that the additional terms in the expression for (A) are not 
generally small, because some denominators in the second sum become extremely small 
and the corrections to the FLRW distance become of order MjL^ or even 1, instead 
of being of order M/L; see [3] for a detailed discussion of these effects. By carefully 
studying these additional terms, we arrive at the conclusion that the perturbative 
corrections to the FLRW distance/redshift relation remain small provided: 

This relation between the mass of the object, M/L, and the compacity of the lattice, 
?7/L, shows that if objects are too compact, the perturbative expansion breaks down, 
as far as the calculation of observables down a past lightcone is concerned, even though 
the perturbative calculations remain a good way of estimating the spacetime geometry 
(i.e. of solving the Einstein field equations). Similar problems were encountered 
in perturbations of an FLRW background in [7]. If this criterion is not satisfied, 
the perturbative calculations cannot be trusted, as second order terms (in a/ M/L) 
become of order 0{\\. one needs to integrate the system of Sachs equations without 
any perturbative expansion, thus retaining the coupling between isotropic expansion 
and shear. This is an ongoing work ^ and preliminary results indicate that when 
Eq. (I7j) is not satisfied, the contribution of the shear modifies significantly the FLRW 
observables: the corrections are usually smaller than the order 1 corrections predicted 
by the perturbative expansion presented here, but they are definitely significant to 
raise the issue of a fitting problem. For example. Fig. [1] shows the percent change 
5r^(z) = 100 X (r^(-2) — r^^^'^(z)) /r^^^'^(z) in the angular distance between a lattice 
with M = IO^^Mq, L = 1 Mpc, ri = O.OIL (lattice of galaxy-like objects) that 
does not satisfy the criterion ([7]) and the equivalent smoothed FLRW model, obtained 
from a complete integration of Sachs equations. We see that the 'divergence' problem 
encountered for such lattices when using perturbative methods is somehow 'cured' by 
solving the full system, even though, differences appear and seem to be systematically 
increasing with the redshift, irrespective of the direction on the sky. 

Table [1] presents the typical values of M, L and rj for some lattices of standard 
astrophysical objects and shows whether the criterion (I7l) is satisfied or not for such 
lattices. 

We see that a lattice of galaxies composed of their disk only marginally fails to 
pass the criterion, whereas when we include the Dark Matter halo, they pass the test 
marginally. A lattice of cluster-like objects passes the criterion easily, but rj/L ~ 1, and 
one can hardly talk of a lattice of separated objects in that case: such objects could not 
really be considered as viriahsed, independent objects as we did in this work. 
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Figure 1. Percent change SrA{z) in tlie angular distance between a lattice with 
M = IO^^Mq, L = 1 Mpc, rj = O.OIL and its FLRW counterpart, obtained from 
a complete integration of Sachs equations. The result is presented for 25 different 
directions on the celestial sphere of the observer located at the centre of one of the 
masses. 



Object 


Rs 

(Mpc) 


L 

(Mpc) 


V 

(Mpc) 


Criterion 
passed 


Neutron 
star 


10-19 


10-6 


10-18 


No 


Galaxy 
(disk) 


10-8 


1 


10-2 


No 

(Marginally) 


Galaxy 
(disk+DM halo) 


10-7 


1 


0.05 


Yes 
(Marginally) 


Galaxy 
cluster 


10-4 


30 


20 


Yes 



Table 1. Some typical lattices and their characteristic parameters. The choices are 
only indicative. The last column answers the question: does such a lattice satisfy the 
criterion ([7])? 

3. Discussion 

We have presented a toy model of the Universe in the form of a regular cubic lattice of 
equal masses of typical size t] whose kinematics is identical on large scales to the FLRW 
model obtained by smoothing the distribution of masses; this model does not exhibits 
any backreaction. We have seen that, despite the fact that a perturbative expansion in 
terms of Rs/L gives a very good approximation for the geometry of space-time, it is 
not suitable for the accurate calculation of observables in the model if the objects are 
too compact. Specifically, we have shown that a perturbative calculation of observables 
can be trusted only if the parameters of the lattice satisfy ([7]). If this is not the case, 
a non perturbative approach is needed to fully take into account the effect of the Weyl 
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curvature sourcing the shear of bundles of null geodesies. Somehow, this was to be 
expected: the perturbative calculation decouple the shear from the isotropic expansion, 
making the observables independent on the Weyl curvature, but we know that in a 
mostly empty Universe (masses very compact), the behaviour of null ray bundles must 
be dominated by the Weyl curvature. The bound (171) gives a quantitative criterion to 
decide what 'too compact' means. This also illustrates the importance of the fitting 
problem: if ([7]) is not satisfied, the FLRW reconstructed by smoothing spatially the 
kinematics of the model differs systematically and sometimes significantly from the 
FLRW model fitting the observations on the past lightcone of observers. Choosing the 
order of magnitude of the parameters of the model to represent something 'realistic' is 
difficult, but it is interesting to note from Tabled] that galaxy-like objects and lattices 
are exactly at the transition between the lattices that pass the criterion ([7]) and those 
which do not. This might be extremely important in the precise characterisation of 
the properties of Dark Energy, and the detailed exploration of the consequences of this 
bound as well as the precise non perturbative estimates in the cases when it is not 
satisfied are the subject of an ongoing work [8]. 
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